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INTRODUCTION 


Non-linear  models  arise  in  many  practical  applications . Models 
based  on  Mlchaelis-Menten  differential  equations  are  commonly  used  In 
biochemistry,  biopharmaceutics,  pharmacokinetics  and  other  related 
disciplines.  These  differential  equations  result  in  general  when  com- 
partmental  models  are  applied.  Linear  differential  equations  usually 
suffice  for  the  study  of  compartmental  models  with  constant  rates  of 
exchange  among  compartments.  However,  there  are  many  situations  espec- 
ially in  the  study  of  drug  distribution  where  nonlinearity  arises 
naturally.  Several  methods  have  been  proposed  in  the  literature  to  fit 
constants  In  such  nonlinear  models.  In  an  earlier  paper,  Rustagi  and 
Singh  (1977)  used  difference  equation  approach  to  fit  one-  and  two- 
compartment  models  under  the  assumptions  of  linear  kinetics. 

In  this  paper,  the  Michaelis-Menten  elimination  scheme  is 
studied  using  the  difference  equation  approach.  The  distribution  of 
the  estimates  of  the  resulting  rate  constants  is  not  easily  obtainable 
in  closed  form.  Using  Monte-Carlo  methods,  the  distribution  of  these 
constants  is  derived  for  the  special  case  of  normal  errors.  Similar 
methods  can  be  utilized  under  different  distributional  assumptions 
for  errors  in  the  model. 

Applications  of  Mlchaelis-Menten  equations  to  pharmacokinetics 
has  recently  been  made  by  Wagner  (1973)  and  Sedaan  and  Wagner  (1974) . 
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They  obtained  estimates  of  the  rate  constants  by  numerical  methods  but 
did  not  provide  the  distribution  of  these  estimates.  It  is  well  known 
that  It  is  not  easy  to  obtain  the  distribution  of  estimates  which  are 
obtained  through  numerical  techniques.  These  distributions  are,  how- 
ever, necessary  for  statistical  analysis,  for  example,  in  obtaining  the 
confidence  Interval  estimates  for  the  rate  constants  or  in  comparing 
two  different  rate  constants,. 


2.  MICHAELI S-MENTEN  MODELS 


In  a study  dealing  with  enzyme  kinetics , Mlchaells  and  Menten 
(1913)  provided  a model  for  the  following  one  compartment  system, 
given  by  Figure  2.1. 


Figure  2.1 


C is  the  concentration  of  certain  substance  in  the  compartment,  K is 

IS 

the  rate  of  output  and  V is  the  maximum  rate  of  reaction.  This  system 
can  be  represented  by  the  differential  equation 


dC 

dt 


V c 

m 


K + C * 


m 


(2.1) 


Consider  the  difference  equation  analog  of  the  equation  (2.1)  assuming 
that  the  process  is  observed  at  times  t^.t^, > • • .t^.  Using  the 
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approximation  for  derivatives  In  the  case  of  unequal  time  intervals,  we 


c(ti+i)~c(ti)  v(ti> 


ci+l  “ ci 


Km  + c<ci>  * 


(2.2) 


Equation  (2.2)  is  simplified  in  the  following  form. 


C(t  )-C(t  ) C(t  )-C(t  ) 

K -L-  + V C(t.)  + C(t,)— - £2 -i-  - 0 

m tm  - t mi  1 t - t 


(2.3) 


The  model  for  estimating  constants  K and  V is  assumed  to  be 

mm 


C(t  )-C(t  ) C(t,..)-C(t.) 

n - " + y -C(t  ) + C(t  ) -i-  - u(t.)  . (2.4) 

Vr'i  * 1 1 ‘i+i  “ ci  1 


We  assume  for  simplicity  that  u(t^)  are  random  errors  with  zero  means 
2 

and  variances  . 

The  weighted  least-squares  solution  for  K and  V are  obtained 

nx  m 


by  minimizing  the  expression! 


S(V  »X  ) - I ■—  [k  C(ti+l)“C(ti)  v f c(t1+1)-c(ti)j 

m m H TT + vmc(c±)  + Z ~ 


1-1  1 I 


1 *!+!-*! 


where  — is  the  weight  for  each  data  point*  By  taking  the  derivatives 
i 

of  equation  (2.5)  with  respect  to  K and  V , we  get  the  normal  equations 

m t& 

given  in  (2.6)  end  (2.7)  which  provide  the  least-squares  estimates  for 

K and  V . 
m m 
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a n-1 
K £ 


_1_  (C(ti-fi:>~C<ti? 


“ i"l  ai  \ 


+ V £ ic(C  i- 

B i-i  ai  1 ti+rt± 


n-1  C(t  ) / C(t...)-Ct.)  2 

+ £ I ■ ■ ■ . — i-  - 0 


1-1  *1 


ti+l“ti 


~ n-1  C(t.)  C(t...)-Ct.)  ^ n-1  (C(t.))w 

r r ...  . aZs  — + tr  r L 

“ l-i  ai  ci+i_ti  ® i-i  ci 


n-1  (COO)' 

+ £ -i — 


1-1 


c(ti+1)-c(ti) 

tl+l"Ci 


(2.6) 


(2.7) 


Example 

Suppose  the  time-concentration  data  of  alcohol  elimination  in 
human  subject  is  given  In  Table  2.1.  The  data  is  fitted  to  (2.4)  by 
using  — — - — as  the  weights,  we  get  K - 2.8021  mM,  and 

(c(tl))2  " 

V - 0.0882  mM/min.  The  ratio  V /K  is  0.0315  min  This  ratio  is 
m rn  m 

close  to  the  theoretical  approximate  first-order  rate  constant  for  drug 
elimination  which  is  known  to  be  0.0513  min~\  Wagner  (1971). 


Table  2.1 


Time  (min) 

5.0 

48.0 

78.0 

105.0 

135.0 

163.0 

Concen.  (mM) 

6.9 

4.1 

2.7 

1.4 

0.5 

0.15 

Note  that  the  distribution  of  the  estimates  K and  V , which 

m m 

depend  non-linearly  on  C(t  ),  are  not  easily  obtainable  for  given  n. 
Simulation  methods  are  used  later  to  obtain  the  properties  of  these 


estimates. 
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Char*ct«*^«  of  the  twc-compartnent  model  wlch 
MlchaeXfs-Haatea  have  ataiM  hy  f- 

■*'  “*8""  «”«>•  8“PP<>«  • in,  i.  f.J.e,*  lnt0 

• ^ a?.te.  fatr.vea.aa1?.  ,Mte.  „ b<  ctpr,sw<d  b?  _ 

“d'1  *«  **«-«-  fr„  the  flr.t  c„. 

18  *»  »•»  2.2  aa^ 

kinetics. 


c2Ct) 


, ^(t) 

IT  [ v 

a m 


Figure  2.2 

Cl(W  “d  C2(t)  “*  th*  «oM«».tlon.  fa  the  tvo  c<»p.rt»eat., 

““  - h21  are  th.  teat-order  tat.  caataata.  .ad  , a...  th.  a™ 

—1.S  a.  before.  The  a.th.»tlt.l  aod.X  far  taaceatt.tfoaa  re..!..  f. 
the  equations  (2.8)  and  (2.9). 


21  "V^0!  + k12C2 


(2.8) 


dt  * k21Cl  - ki2C2  * 


Suppose  an  experimenter  can  observe  the 
partnent.  From  (2.8)  we  have 


concentration  in  only  one  com- 


UVM  V 

k12C2  “ 7T  + k2ici  + F-fc^i 

B 1 


(2.10) 


a..-*--—...,.,-  a!.-\  '*,t-SKv.  , i .tS 


M 


l 
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L 


Substituting  (2.10)  into  (2.9) , we  get 


V 

m 

K + C 


a 


1 


6 


(2.11) 


Differentiating  equation  (2.8)  once  more,  we  have  a second-order  dif- 
ferential equation  in  terms  of  C^, 

2 2 dCl 

+l(W  ft21-k12)-W-dr+k12VK.+Cl)Cl'°-  (2-12> 


Suppose  now  that  the  concentration  in  a given  compartment  Is 

observed  at  equal  intervals  of  time.  For  notatlonal  conveni-nce  we 

dC 

drop  the  subscript  in  and  replace  the  derivative  — by  C(t+1)  - C(t) 
.2 

and  — s-  by  C(t+2)  - 2(t+l)  + C(t)  to  obtain  the  difference  equation 
dt 

analog  to  equation  (2.12).  We  have 

(C(M-2)-  2C(t+l)+C(t)][K  +C(t)]2+  [C(t+1)  -C(t)]{[K  +C(t)]2 

m is 

<k2rk12)-Vm}  + C<t)[Km+C(t)lk12Vm  * 0 ' (2*13> 

The  model  for  estimating  the  rate  constants  is 

[C(t+2)  - 2C(t+l)  +C(t)][K  + C(t)  ] 2 + [C(t+1)  -C(t)]{[K  +C(t)]2 

m n 

(k-.-k. ,)  - V K } + C(t) (K  + C(t) ]k1 -V  - u(t)  (2.14) 
LX  XL  m m m XL  m 

where  u(t)  is  the  error  term  with  zero  mean  and  appropriate  variance. 

Notice  that  the  equation  (2. 14)  is  a polynomial  with  degree 

three  in  parameters  K , V , lc,-,  and  k_,  and  results  in  a complicated 

m ifl  xl  lx 
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estimation  procedure.  With  certain  additional  assumptions,  such  as 
those  of  some  parameters  known,  nonlinear  least-squares  method  may  be 
utilized.  The  uniqueness  of  the  resulting  estimates  is  still  question- 
able. 


3.  ESTIMATION  OP  BATE  CONSTANTS 

Several  methods  of  estimation  of  rate  constants  in  pharmaco- 
kinetics models  are  in  common  use.  Numerical  procedures  leading  to 
techniques  using  quasilinearization  are  due  to  Bellman  and  Kalaba 
(1965) . Also  computer  programs  such  as  NONLXN,  by  Ketzler  (1969)  which 
is  based  on  the  modification  of  Gauss-Newton  procedure  are  commonly 
used.  We  use  the  difference-equation  approach  discussed  earlier  for 
obtaining  the  estimates. 

Using  Runge-Kutta  method,  which  is  a part  of  NONI.IN  program, 

data  are  generated  from  the  Kichaells-Kenten  model  (2.1)  with  a given 

set  of  values  of  constants  K end  V . A sample  of  25  observations  is 

mm 

generated  and  each  value  is  increased  by  a standard  normal  random  vari- 
ate. The  estimates  of  V and  K are  then  made  using  the  difference 

mm 

equations  model  of  the  Michaelis-Menten  equations  as  in  equation  (2.4), 

The  resulting  estimates  are  given  in  Table  I for  one  hundred  such  samples. 

A A 

The  frequency  distributions  of  K and  V are  given  in  Tables  II 

m in 

A A 

and  III.  Means  and  standard  deviations  for  K and  V are 

tn  m 
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To  test  the  hypotheses  that  the  sample  of  one  hundred  K and 

A 

A 

V arise  from  normal  distributions,  classical  chi-squared  test  of 

ft 

goodnesa-of-f it  and  Lilliefors  test  are  used,  Conover  (1971).  The 

A A 

chi-squared  values  for  the  samples  K and  V are 

10  A 

2 

Xg  m 6.454  with  5 degrees  of  freedom, 
a 
2 

Xy  “ 2.516  with  4 degrees  of  freedom, 
a 

which  are  not  significant  at  0.05  level. 

For  Lilliefors  ncn-parometric  test,  the  test  statistic  used  is 

T - Sup  |F(x)  - S(x)|  , 
x 

where  F(x)  is  the  standard  normal  distribution  function  and  S(x)  is 
the  empirical  distribution  function.  The  values  have  been  standardized 
using  the  means  and  standard  deviations  obtained  above.  Absolute  values 
of  (F(x)-S(x))  for  both  samples  are  recorded  in  Table  Z.  The  maximum 

A 

value  for  K occurs  at  x - 0.23,  where  S(x)  equals  0.53,  F(x)  equals 
111 

A 

0.59,  and  T is  0.06.  The  maximum  value  for  V occurs  at  x - -0.56, 

m 

where  S(x)  equals  0.35,  F(x)  is  0.29  so  that  T equals  0.06.  The  maxi- 
mum value  0.06  also  occurs  at  other  points,  but  at  no  point  does  the 
absolute  difference  of  S(x)  and  F(x)  exceed  0.06. 

The  Lilliefors  test  calls  for  rejection  of  our  hypotheses  at 
a » 0.2C  if  T exceeds  its  0.80  quantile.  The  critical  region  is 
obtained  from  tables  in  Conover  (1971).  The  hypotheses  are  accepted, 

A A 

and  therefore  we  conclude  that  the  estimates  K and  V are  normally 

1ft  Ill 


distributed. 
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0.0 

-1.24 

0.0 

-1.24 

0.0 

-1.24 

0.0 

-1.24 

0.0 

-1.24 

0.0 

-1.24 

0.0 

-1.24 

0.0 

-1.24 

0.0005 

-1.23 

0.0027 

-1.18 

0.0045 

-1.14 

0.0055 

-1.12 

0.0058 

-1.11 

0.0078 

-1.07 

0.0100 

-1.02 

0.010 

-1.02 

0.0114 

-0.99 

0.0123 

-0.97 

0.0204 

-0.79 

0.0206 

-0.78 

0.0221 

-0.75 

0.0221 

-0.75 

0.0225 

-0.74 

0.0243 

-0.70 

0.0262 

-0.66 

0.0295 

-0.59 

0.0299 

-0.58 

0.0305 

-0.56 

0.0309 

-0.56 

0.0318 

-0.54 

0.0325 

-0.52 

0.0332 

-0.50 

0.0344 

-0.48 

0.0357 

-0.45 

0.0362 

-0.44 

0.0377 

-0.40 

0.0379 

-0.40 

0.0394 

-0.37 

0.0426 

-0.30 

0.0429 

-0.29 

0.0429 

-0.29 

0.0439 

-0.27 

0.0449 

-0.25 

0.0497 

-0.14 

0.0499 

-0.13 

0.0501 

-0.13 

0.0503 

-0.13 

0.0531 

-0.06 

.02 

0.0415 

.02 

0.1407 

.02 

0.1508 

.02 

0.1567 

.02 

0.1611 

.02 

0.1665 

.02 

0.1797 

.02 

0.1829 

.01 

0.1864 

.01 

0.1875 

.01 

0.1882 

.01 

0.1883 

.003 

0.1892 

.002 

0.1895 

.01 

0.1915 

.01 

0.1938 

.01 

0.1945 

.02 

0.2020 

.02 

0.2037 

.01 

0.2039 

.006 

0.2044 

.006 

0.2056 

.01 

0.2061 

.002 

0.2072 

.004 

0.2095 

.01 

0.2101 

.01 

0.2102 

.01 

0.2107 

.01 

0.2108 

.01 

0.2130 

.01 

0.2130 

,02 

0.2130 

.02 

0.2131 

.02 

0.2135 

.02 

0.2144 

.03 

0.2162 

.03 

0.2205 

.03 

0.2210 

.01 

0.2217 

.03 

0.2232 

,03 

0.2264 

,03 

0.2329 

,03 

0.2330 

004 

0.2330 

03 

0.2350 

03 

0.2354 

03 

0.2360 

01 

0.2360 

-4.05 

.01 

-2.05 

.002 

-1.84 

.002 

-1.72 

.002 

-1.63 

.001 

-1.53 

.003 

-1.26 

.03 

-1.19 

.04 

-1.12 

.04 

-1.10 

.04 

-1.09 

.03 

-1.08 

.02 

-1.07 

.01 

-1.06 

.004 

-1.02 

.003 

-0.97 

.006 

-0.96 

.002 

-0.81 

.03 

-0.77 

.02 

-0.77 

.02 

-0.76 

.01 

-0.74 

.009 

-0.73 

.003 

-0.70 

.0  02 

-0.66 

.005 

-0.64 

.01 

-0.64 

.01 

-0.63 

.03 

-0.63 

.03 

-0.59 

.04 

-0.59 

.04 

-0.59 

.04 

-0.58 

.06 

-0.58 

.06 

-0.56 

.06 

-0.52 

.05 

-0.43 

.04 

-0.42 

.04 

-0.41 

.05 

-0.38 

.05 

-0.32 

.04 

-0.18 

.01 

-0.18 

.01 

-0.18 

.01 

-0.14 

.006 

-0,13 

.01 

-0.12 

.040 

-0.12 

.04 
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TABLE  I (continued) 


K (K  -.056)/. 045  S(x)-?(x)  V (V  -.242)/. 0495  S(X)-F(X) 

la  in  m . 


0.0558 

-0.00 

.01 

0.0567 

0.02 

.008 

0.0583 

0.05 

.009 

0.0610 

0.11 

.02 

0.0662 

0.23 

.06 

0.0667 

0.24 

.02 

0.0668 

0.24 

.02 

0.0668 

0.24 

.02 

0.0670 

0.24 

.02 

0.0674 

0.25 

.02 

0.0678 

0.26 

.01 

0.0689 

0.29 

.01 

0.0699 

C.31 

.01 

0.0740 

0.40 

.03 

0.0758 

0.44 

.04 

0.0764 

0.45 

.03 

0.0769 

0.46 

.02 

0.0805 

0.54 

.03 

0.0806 

0.54 

.03 

0.0824 

0.58 

.03 

0.0827 

0.59 

.03 

0.0841 

0.62 

.03 

0.0844 

0.63 

.02 

0.0853 

0.65 

.02 

0.0880 

0.71 

.03 

0.0885 

0.72 

.02 

0.0895 

0.74 

.02 

0.0906 

0.77 

.01 

0.0906 

0.77 

.01 

0.0909 

0.77 

.01 

0.0925 

0.81 

.001 

0.0946 

0.85 

.002 

0.0952 

0.87 

.01 

0.0984 

0.94 

.006 

0.0994 

0.96 

.001 

0.0998 

0.97 

.01 

0.1001 

0.98 

.02 

0.1028 

1.04 

.01 

0.1043 

1.07 

.02 

0.1090 

1.17 

.01 

0.1137 

1.28 

.01 

0.1138 

1.28 

.01 

0.1143 

1.29 

.01 

0.1169 

1.35 

.01 

0.1178 

1.37 

.02 

0.2360 

-0.12 

.04 

0.2375 

-0.09 

.04 

0.2378 

-0.08 

.04 

0.2415 

-0.01 

.02 

0.2435 

0.03 

.02 

0.2440 

0,04 

.03 

0.2445 

0.05 

.03 

0.2450 

0.06 

.04 

0.2487 

0.14 

.02 

0.2539 

0.24 

.01 

0.2573 

0.31 

.03 

0.2585 

0.33 

.03 

0.2595 

0.35 

.03 

0.2601 

0.37 

.02 

0.2625 

0,41 

.03 

0.2673 

0.51 

.05 

0.2675 

0.52 

.05 

0.2707 

0.58 

.06 

0.2710 

0.59 

.05 

0.2741 

0.65 

.06 

0,2747 

0.66 

.05 

0.2759 

0.68 

.050 

0.2761 

0.69 

.05 

0.2774 

0.72 

.04 

0.2784 

0.74 

.03 

0.2788 

0.74 

.03 

0.2809 

0.79 

.02 

0.2810 

0.79 

.02 

0.2812 

0.79 

.02 

0.2818 

0.80 

.008 

0.2824 

0.82 

.004 

0.2852 

0.87 

.008 

0.2873 

0.92 

.01 

0.2884 

0.94 

.006 

0.2892 

0.95 

.01 

0.2892 

0.95 

.01 

0.2914 

1.00 

.009 

0.2923 

1.02 

.01 

0.2974 

1.12 

.003 

0 . 3006 

1.18 

.001 

0.3032 

1.24 

.003 

0.3052 

1.28 

.001 

0.3057 

1.29 

.02 

0.3058 

1.29 

.01 

0.3094 

1.36 

.02 
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TABLE  I (continued) 


K 

m 

(K  -.056)/  .045 

SB 

S(x)-F(x)  V 

So 

(V  -.242)/. 0495  S(X)-F(X) 

Q 

0.1201 

1.42 

.02 

0.3114 

1.40 

.02 

0.1206 

1.43 

.03 

0.3161 

1.50 

.02 

0.1222 

1.46 

.04 

0.3233 

1.64 

.01 

0.1232 

1.49 

.04 

0.3248 

1.67 

.02 

0.1303 

1.64 

.04 

0.3250 

1.68 

.03 

0.1311 

*.66 

.04 

0.3283 

1.74 

.04 

0.1322 

1.69 

.05 

0.3333 

1.84 

.04 

13 


TABLE  II 

FREQUENCY  DISTRIBUTION  OF  K 

in 


A 

K 

m 

Frequencies 

Relative 

Frequencies 

Cumulative 

Relative 

Frequencies 

0 

8 

0.08 

0.08 

0.00-0.02 

10 

0.10 

0.18 

0.02-0.04 

20 

0.20 

0.38 

0.04-0.06 

13 

0.13 

0.51 

0.06-0.08 

14 

0.14 

0.65 

0.08-0.10 

19 

0.19 

0.84 

0.10-0.12 

9 

0.09 

0.93 

0.12-0.14 

7 

0.07 

1.0 

Total 


100 


1.0 
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I 


i 

i 


TABLE  III 

FREQUENCY  DISTRIBUTION  OF 


A 

I 

n 

Frequencies 

Relative 

Frequencies 

Cumulative 

Relative 

Frequencies 

0.00-0  17 

6 

0.06 

0.06 

0.17-0.20 

11 

0.11 

0.17 

0.20-0.23 

24 

0.24 

0.41 

0.23-0.26 

20 

0.20 

0.61 

0.26-0.29 

23 

0.23 

0.84 

0.29-0.32 

11 

0.11 

0.95 

0.32-0.35 

5 

0.05 

1.0 

Total 

100 

1.0 
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